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The  present  report  summarizes  the  research  performed  under  Contract  DAAG29-84- 
K-0191  on  various  problems  related  to  visibility  probabilities,  when  the  fields  of  vision 
are  obscurred  by  random  elements.  Several  Technical  Reports  have  been  published.  The 
paper  on  the  three-dimensional  problem  was  published  in  the  Naval  Research  Logistics 
Quarterly  (1988).  The  paper  on  the  multi-observer  multi-target  visibility  probabilities 
is  under  revision  for  the  Naval  Research  Logistics  Quarterly.  The  Principal  Investigator 
has  also  written  a  monograph  on  the  subject  based  on  a  workshop  manual.  The  research 
problems,  the  Table  of  Contents  for  the  monograph  and  the  introductory  sections  of  the 
Technical  Reports  follow  hereby. 
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THE  VISIBILITY  OF  STATIONARY  AND  MOVING  TARGETS  IN  THE 
PLANE  SUBJECT  TO  A  POISSON  FIELD  OF  SHADOWING  ELEMENTS 


M  N  A  DIN.*  Technton  —  Israel  Institute  of  Technology 
S  ZaCKS.’*  State  University  of  \Vh  )ork  at  Binghamton 


abstract 

A  methodology  for  an  analytical  derivation  of  visibility  probabilities  of  n 
stationary  target  points  in  the  plane  is  developed  for  the  case  when  shadows  are 
cast  b>  a  Poison  random  field  of  obscuring  elements  In  addition,  formulae  for 
the  moments  of  a  measure  of  the  total  proportional  visibility  along  a  star-shaped 
curve  are  given. 

VISIBILITY  PROBABILITY.  MEASURE  OF  PROPORTIONAL  VISIBILITY.  POISSON 
RANDOM  SHADOWING  PROCESS 


0.  Introduction 

.present  paper  discusses  problems  of  visibility  of  targets  through  randefm 
fields  of  Dbscurmg  elements  (trees,  bushes,  clouds,  etc.).  The  problems  dealt  with 
are  of  a  stochastic  nature.  The  exact  number,  location  and  dimensions  of  the 
obscuring  elements  are  unknown.  A  probability  model  is  formulated  concerning 
these  variables.  Given  such  a  probability  model,  it  is  required  to  determine 
probabilities  of  certain  events  and  distributions  of  certain  random  variables.  To 
illustrate  some  of  the  problems  that  can  be  solved  bv  the  methodology 
developed,  consider  the  following  examples. 


r 

/ 


Example  1.  \'isibility  of  stationary  targets.  An  observer  is  placed  at  a  given 
location  in  a  forest,  in  order  to  detect  specified  targets  (vehicles,  animals,  etc  ). 
Due  to  the  random  location  of  the  trees  it  is  important  to  determine  the 
probabilities  that  individual  targets  are  observed  and  the  distribution  of  the 
number  of  targets  observed.  For  this  purpose  one  has  to  determine  the 
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Visibility  Probabilities  on  Line  Segments  in 
Three-Dimensional  Spaces  Subject  to  Random 
Poisson  Fields  of  Obscuring  Spheres* 


M.  Yadin  and  S.  Zacks 

Technion,  Israel  Institute  of  Technology  and  State  University  of  New  York  at 

Binghamton 


The  methodology  of  determining  simultaneous  visibility  probabilities  and  points  on  line 
segments  is  extended  to  problems  of  three-dimensional  spaces,  with  Poisson  random  fields 
of  obscuring  spheres.  Required  functions  are  derived  analytically  and  a  numerical  example 
is  given  for  a  special  case  of  a  standard  Poisson  field,  with  uniform  distribution  of  sphere 
diameters. 


1.  INTRODUCTION 

The  present  article  studies  a  three-dimensional  shadowing  problem,  and  provides 
an  algorithm  for  the  computation  of  visibility  probabilities,  for  cases  where  the  target 
is  a  straight  line  segment,  or  points  on  a  line  segment,  which  is  parallel  to  a  layer  of 
shadowing  objects.  The  results  given  here  can  be  generalized  to  target  curves  in  three- 
dimensional  space,  provided  they  lie  in  a  plane  containing  the  observation  point  O. 
As  we  show  in  Section  5,  the  theory  developed  here  can  be  applied,  for  example,  to 
determine  visibility  probabilities  for  a  pilot  approaching  a  runway  under  cloudy  con¬ 
ditions.  Many  other  applications  are  possible  for  military  and  civil-type  problems. 

Consider  a  problem  of  visibility  obstruction,  when  an  observer  (source  of  light)  is 
located  at  the  origin  O  and  the  obscuring  (shadowing)  elements  are  randomly  dispersed 
between  the  origin  and  the  target  curve.  A  straight  line  segment  is  specified  in  a 
three-dimensional  Euclidean  space.  The  obscuring  elements  are  modeled  as  spheres 
whose  centers  are  randomly  dispersed  within  a  layer  bounded  between  two  planes 
parallel  to  c€,  between  the  origin  and  the  target  curve.  It  is  further  assumed  that  the 
centers  of  spheres  are  dispersed  according  to  a  homogeneous  Poisson  process  in  the 
layer.  The  assumption  about  the  homogeneity  of  the  Poisson  process  is  not  crucial 
and  can  be  relaxed.  The  radii  of  spheres  are  i.i.d.  random  variables  having  a  distribution 
independent  of  the  locations  of  the  centers.  The  problem  is  to  determine  simultaneous 
visibility  probabilities  for  n  arbitrarily  specified  points  on  '€.  Moments  of  integrated 
measures  of  visibility  on  ^  can  be  determined,  as  shown  in  (2),  by  integrating  these 
visibility  probabilities  functions. 

The  methodology  of  the  present  study  is  to  reduce  the  three-dimensional  shadowing 
problem  to  a  two-dimensional  one,  on  the  plane  jW  containing  the  line  segment  %  and 
the  origin  O.  The  treatment  of  the  shadowing  problem  on  M  can  follow  the  methodology 
developed  in  our  previous  studies  (4,5).  There  are,  however,  several  complications. 
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ABSTRACT 

A  shadowing  process  in  the  plane  is  studied,  with  respect  to  a  single  source  of  light 
and  a  Poisson  field  of  shadowing  elements,  located  between  the  source  of  light  and  a 
target  curve.  N  subsegments  of  length  6  are  considered  within  a  specified  portion  of  the 
target  curve.  Random  measures  of  visibility  and  random  shadow  weights  on  the  target 
curves  are  approximated  by  corresponding  discrete  random  variables,  defined  for  the  N 
subsegments.  Algorithms  are  provided  for  the  computation  of  the  distribution  functions 
of  the  approximating  discrete  random  variables. 


Key  words:  Shadowing  processes,  Semi-Markov  processes,  discrete  approximation,  Poisson 
fields,  measure  of  visibility. 
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1.  INTRODUCTION. 

Consider  a  source  of  light  and  a  linear  path  in  a  plane.  Random  objects  are  scattered 
between  the  source  of  light  and  the  path ,  casting  shadows  on  portions  of  it. 

The  Shadovnng  Model  was  originally  introduced  by  Chemoff  and  Dally  [l].  In  their  paper 
they  assumed  that  the  shadowing  objects  are  disks  having  centers  which  are  distributed 
in  a  rectangle  between  the  source  of  light  and  the  path,  according  to  a  two  dimensional 
Poisson  Process,  and  their  diameters  are  represented  by  i.i.d.  random  variables.  The 
paper  applies  queueing  theoretical  techniques  to  derive  the  conditional  distributions  of  the 
lengths  of  the  shadowed  and  the  unshadowed  segments  of  the  path. 

The  shadowed  and  unshadowed  segments  of  the  path  are  analogous  to  busy  and  idle 
periods  in  certain  queueing  processes  with  Markov  arrivals.  Like  in  these  processes,  the 
derivation  of  the  distribution  of  the  lengths  of  unshadowed  segments  (which  are  similar  to 
idle  periods)  is  straightforward.  Chemoff’s  and  Daily’s  paper  concentrates  on  the  distri¬ 
butions  of  shadowed  segments,  and  provide  some  integral  equations  for  these  distributions. 
The  solution  of  these  equations  requires  generally  lengthy  numerical  procedures. 

Yadin  and  Zacks  (2,  3]  introduced  the  concept  of  Random  Measures  of  Visibility,  which 
are  defined  as  the  total  length  of  the  unshadowed  subsegments  of  a  specified  segment  of 
the  path,  or  the  proportion  of  the  segment  which  is  in  the  light.  In  these  papers  the 
distributions  of  Measures  of  Visibility  were  approximated  by  a  mixture  of  a  two  points 
distribution  on  {0, 1}  and  a  beta  distribution  on  (0, 1) .  These  mixed  beta  distributions 
were  obtained  by  equating  the  first  three  moments  of  the  distributions  of  Measures  of 
Visibility  to  the  first  three  moments  of  a  mixed-beta  distribution.  In  [2]  the  path  was 
assumed  to  be  a  segment  of  a  circle,  centered  at  the  origin  at  which  the  source  of  light  was 
located.  The  centers  of  the  disks  were  assumed  to  be  distributed  within  an  annular  region 
in  the  circle.  The  results  were  extended  in  (3j  and  [4],  which  considered  more  general  paths 
and  fields  of  obscuring  objects.  The  methodology  was  extended  to  a  three  dimensional 
model  in  [5j.  Certain  applications  were  discussed  in  (7). 
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The  main  objective  of  the  present  paper  is  to  replace  the  mixed  beta  distributions  by 
discrete  distributions  which  approximate  the  required  distributions.  More  specifically,  we 
approximate  the  distribution  of  a  measure  of  visibility,  by  the  distribution  of  a  discrete 
random  variable,  being  the  number  of  the  proportion  of  the  unshadowed  subsegments, 
which  are  properly  located  on  some  segment  of  the  path. 

The  shadowing  model  and  the  approximated  measures  of  visibility  and  lengths  of  shad¬ 
owed  segments  are  formally  defined  in  section  2.  The  approximating  distributions  can  be 
given  in  terms  of  probabilities  of  elementary  events,  which  are  the  events  signifying  which 
of  the  N  subsegments  are  completely  unshadowed  among  the  set  of  all  N  subsegments. 
These  probabilities  can  be  obtained  from  probabilities  of  simultaneous  visibility,  namely 
the  probabilities  that  all  subsegments  in  a  specified  sets  are  unshadowed.  Simultaneous 
visibility  probabilities  were  determined  for  various  cases  in  [4].  Appendix  A  presents  the 
methodology  needed  for  the  derivations  of  probabilities  of  simultaneous  visibility. 

The  accuracy  of  the  discrete  approximation  to  the  distributions  considered  depends  on 
the  number  of  subsegments.  Assuming  a  set  of  N  subsegments,  one  has  to  compute  2N 
probabilities  of  elementary  events  from  the  same  number  of  probabilities  of  simultaneous 
visibility. 

Any  algorithm  which  is  based  on  the  probabilities  of  simultaneous  visibility  is  necessarily 
non  polynomial,  thus  is  limited  to  a  relatively  small  number  of  points,  and  may  therefore 
yield  insufficient  accuracy.  Fortunately,  the  distributions  of  the  approximated  measures  of 
visibility  can  be  obtained  by  a  polynomial  algorithm  which  is  discussed  in  section  4. 

The  results  of  a  numerical  example  are  presented  and  discussed  in  section  5.  These 
results  illustrate  the  advantages  and  limitations  of  the  algorithms  developed  in  the  present 
paper  with  respect  to  accuracy,  computer  time  and  memory  requirements. 
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1.  Introduction. 

the  present  paper  we  consider  problems  associated  with  the 
visibility  of  targets  from  several  observation  points,  subject  I 

to  interference  by  a  random  field  of  obscuring  elements.  This 
study  is  an  extension  of  previous  ones  [1,  2,  3,  4],  which  deal 
with  visibility  problems  associated  with  a  single  observation 
point.  We  consider  a  random  field  of  obscuring  elements,  which 
are  centered  in  a  region  located  between  the  observation  points 
and  the  targets.  The  targets  are  stationary  points  in  the  plane. 

Dynamic  targets  will  be  considered  in  another  paper.  The  objec¬ 
tive  of  the  present  study  is  to  develop  efficient  algorithms  for 
the  computation  of  visibility  probabilities.  For  example,  if 
there  are  three  observation  points  and  three  targets,  there  are 
nine  possible  lines  of  sight,  connecting  observation  points  with 
targets.  These  lines  of  sight  might  be  intersected  by  obscuring 
elements.  We  say  that  a  target  T  is  observable  from  0  if  the 
corresponding  line  of  sight  is  not  intersected.  An  elementary 
event  is  a  statement  which  specifies  which  lines  of  sight  are 
intersected  and  which  are  not.  In  the  above  example  there  are 

9 

2  «  512  such  elementary  events.  The  problem  is  to  determine 

the  probabilities  of  elementary  events,  as  functions  of  properties 
of  the  random  field  of  obscuring  elements  -  its  geometrical 
structure  and  its  stochastic  characteristics. 

In  Section  2  we  present  the  shadowing  model,  the  characteristics 
of  the  field,  and  define  a  few  random  variables  of  interest.  For 
example,  the  number  of  targets,  N^,  observable  simultaneously  from 
the  i-th  observation  point;  or  the  number  of  observation  points 
Mj  from  which  the  j-th  target  can  be  simultaneously  observed;  etc. 


3 


Distributions  related  to  these  variables  are  discussed.  In  Section 
3  we  define  simultaneous  visibility  probabilities  of  sets  of  obser¬ 
vation  points/targets;  and  the  probabilities  of  elementary  events. 

We  show  how  vectors  of  elementary  probabilities  can  be  obtained 
from  vectors  of  simultaneous  visibility  probabilities  by  linear 
transformations , using  Hadamard  matrices  of  order  2n.  In  Section  4 
we  develop  the  algorithm  for  the  computation  of  simultaneous  visibility 
probabilities.  Section  5  presents  a  geometrical  development  of 
formulae  required  for  the  determination  of  certain  functions 
defined  in  Section  4.  Finally,  in  Section  6  we  provide  a  numerical 
example,  related  to  three  observation  points  and  three  targets.  In 
this  example  we  illustrate  the  application  of  the  procedures 
developed  in  the  present  study  to  answer  various  questions  of 
interest.  We  tabulate  also  the  joint  probability  distribution 
function  (P.D.F.)  of  (N^,N2»N3),  and  the  three  joint  marginals  of 
pairs  of  random  variables.  These  joint  distributions  reflect  the 
type  of  stochastic  dependence  between  these  random  variables,  and 
show  that  the  regular  correlations  are  very  low.  This  result  is, 
however,  a  function  of  the  particular  location  of  the  observation 
points  and  of  the  targets.  By  changing  the  location  coordinates 
of  these  points  the  correlations  in  the  above  distributions  may  i 

increase  or  decrease.  Different  indexes  of  stochastic  dependence 
can  be  developed  for  different  problems  under  consideration.  The 
simplest  index  of  dependence  can  be  based  on  the  simultaneous  \ 

visibility  probability  of  two  or  more  targets  from  a  given  observa¬ 
tion  point.  Another  index  of  dependence,  could  be  the  correlation 
between  and  N ^  (i#j),  or  and  M ^  (i^j).  The  algorithms  ! 

developed  in  the  present  study  can  serve  as  basic  tools  for  the 
computation  of  such  indexes. 
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1.  Introduction. 

The  present  paper  studies  a  three-dimensional  shadowing  problem,  and  provides  an  algo¬ 
rithm  for  the  computation  of  visibility  probabilities,  for  cases  where  the  target  is  a  straight 
line  segment,  or  points  on  a  line  segment,  which  is  parallel  to  a  layer  of  shadowing  objects. 
The  results  given  here  can  be  generalized  to  target  curves  in  three-dimensional  space,  pro¬ 
vided  they  lie  in  a  plane  containing  the  observation  point  O .  As  we  show  in  Section  5, 
the  theory  developed  here  can  be  applied,  for  example,  to  determine  visibility  probabilities 
for  a  pilot  approaching  a  runway  under  cloudy  conditions.  Many  other  applications  are 
possible  for  military  and  civil  type  problems. 

Consider  a  problem  of  visibility  obstruction,  when  an  observer  (source  of  light)  is  lo¬ 
cated  at  the  origin,  O ,  and  the  obscuring  (shadowing)  elements  are  randomly  dispersed 
between  the  origin  and  the  target  curve.  A  straight  line  segment,  C,  is  specified  in  a 
three-dimensional  Euclidean  space.  The  obscuring  elements  are  modeled  as  spheres  whose 
centers  are  randomly  dispersed  within  a  layer  bounded  between  two  planes  parallel  to  C , 
between  the  origin  and  the  target  curve.  It  is  further  assumed  that  the  centers  of  spheres 
axe  dispersed  according  to  a  homogeneous  Poisson  process  in  the  layer.  The  assumption 
about  the  homogeneity  of  the  Poisson  process  is  not  crucial  and  can  be  relaxed.  The  radii 
of  spheres  are  i.i.d.  random  variables  having  a  distribution  independent  of  the  locations 
of  the  centers.  The  problem  is  to  determine  simultaneous  visibility  probabilities  for  n 
arbitrarily  specified  points  on  C .  Moments  of  integrated  measures  of  visibility  on  C  can 
be  determined,  as  shown  in  [2],  by  integrating  these  visibility  probabilities  functions. 

The  methodology  of  the  present  study  is  to  reduce  the  three-dimensional  shadowing 
problem  to  a  two-dimensional  one,  on  the  plane  X  containing  the  line  segment  C  and  the 
origin  O.  The  treatment  of  the  shadowing  problem  on  M  can  follow  the  methodology 
developed  ia  our  previous  studies  [3,  4j.  There  are,  however,  several  complications.  Even 
if  the  Poisson  field  of  sphere  centers  within  the  layer  is  homogeneous,  the  projection  of  the 
field  into  M  generally  yields  a  non- homogeneous  field  of  disks.  The  radii  of  disks  on  M 


depend  on  the  location  of  the  corresponding  centers  in  the  layer,  the  length  of  the  sphere 
radii  and  their  distances  for  X.  Suppose  that  the  distributions  of  radii  of  spheres  are 
supported  by  a  finite  interval  (0,  b) ,  where  b  is  sufficiently  small,  so  that  neither  the  origin 
O  nor  the  target  C  can  be  covered  by  a  random  sphere.  Obviously,  spheres  having  centers 
with  distance  from  X  greater  than  b  can  be  ignored. 

In  order  to  avoid  computational  difficulties  created  by  the  non-homogeneity  of  the  field 
on  X,  we  partition  the  layer  into  a  large  number  of  thin  sublayers.  The  projection  of 
the  random  fields  restricted  to  these  infinitesimal  sublayers  on  X  are  strips  of  almost 
homogeneous  sub-fields. 

Given  n  specified  points  on  C ,  consider  the  rays  originating  at  O  and  passing  through 
these  points.  The  n  points  are  simultaneously  visibile  (in  the  light)  if  none  of  the  random 
disks  on  X  intersects  any  one  of  these  rays.  Thus,  if  2?(n)  designates  the  expected  number 
of  disks  on  X ,  which  intersect  at  least  one  of  the  specified  rays,  then  the  probability  of 
simultaneous  visibility  is  exp  (~E(n)) .  The  main  contribution  of  the  present  paper  is  in 
developing  the  formulae  for  computing  E(n) .  For  this  purpose  we  present  in  Section  2  the 
three  dimensional  model  and  the  basic  geometry.  In  Section  3  we  discuss  the  reduction  to 
a  two-dimensional  shadowing  process.  In  Section  4  we  present  the  formula  for  computing 
simultaneous  visibility  probabilities.  Analytical  derivation  of  some  of  the  functions  speci¬ 
fied  in  Section  4  is  given  in  the  Appendix.  In  Section  5  we  provide  an  explicit  solution  for 
the  case  of  i.i.d.  sphere  radii  having  a  uniform  distribution  on  (0, 6) . 

In  order  to  assist  the  reader  in  implementing  the  methodology  of  the  present  paper,  we 
highlight  below  the  main  steps  required  for  the  computation  of  visibility  probabilities.  A 
FORTRAN  program  according  to  which  the  example  of  Section  5  has  been  computed  can 
be  furnished  upon  request. 

STEP  X:  Given  O  and  C,  specify  the  parallel  planes  ( U *)  and  (W*)  which  constitute 
the  boundaries  to  the  layer  of  random  spheres;  the  intensity,  A  of  the  Poisson  field  in  the 
layer  and  the  distribution  function  F(-)  of  radii  of  spheres. 
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STEP  2:  Determine  the  plane  C *  containing  C,  which  is  parallel  to  (W*) ;  and  the 
distances  u*,  w* ,  r*  from  O  of  these  planes. 

STEP  3:  Determine  the  projections  O*  and  O"  of  O  on  C*  and  on  C,  respectively. 
Let  d*  be  the  distance  of  O"  from  O* .  The  inclination  angle  of  X  is  <p  =  arctan(d*/r*). 
STEP  4:  Determine  =  6*tan(^)  =  b  ■  d*/r* ; 

STEP  5:  Prepare  a  subroutine  to  compute  the  function  K(s,t)  according  to  formulae 
(4.6),  (4.7)  and  (A. 13).  [The  coordinates  s  are  measured  along  C  relative  to  O" .] 

STEP  6:  Compute  simultaneous  visibility  probabilities,  according  to  (4.1)-(4.3)  and 
(4.9)-(4.10). 

2.  The  Three  Dimensional  Shadowing  Model. 

Consider  a  source  of  light  located  at  a  point  O  (the  origin),  and  an  arbitrary  straight  line, 
C ,  in  a  three-dimensional  Euclidean  space.  A  layer  of  random  spheres  is  located  between 
two  planes  U* ,  W* ,  which  are  parallel  to  the  line  C.  Let  C*  be  a  plane  containing  C, 
parallel  to  U*  (see  Figure  l).  The  random  spheres,  whose  centers  are  located  between  U* 
and  W*  cast  shadows  on  C.  As  stated  in  the  Introduction,  the  objective  is  to  develop 
formulae  for  the  computation  of  the  probabilities  of  simultaneous  visibility  of  any  n  points 
on  C  (n  >  1).  These  probabilities  are  then  applied  to  determine  the  moments  of  total 
visibility  measures  on  C,  and  to  approximate  its  distribution.  In  the  present  section  we 
discuss  the  basic  stochastic  and  geometric  structure  of  the  field. 

Let  u* ,  ti>*  and  r*  be  the  distances  of  U* ,  Ml*  and  C*  from  O,  respectively.  It  is 
assumed  that 

(2.1)  0  <  b  <  u*  <  w*  <  r*  —  6, 

where  6  is  the  maximal  radius  of  a  sphere.  Let  M  be  a  plane  passing  through  O  and  C . 
Let  U  and  W  be  straight  lines  at  which  M  interesects  U*  and  W* .  The  three  parallel 
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